Focusing of high-energy particles 
in the electrostatic field of a homogeneously charged sphere 
and the effective momentum approximation 



^ i Andreas Aste, Dirk Trautmann 

I Department of Physics and Astronomy, University of Basel, 

CN ■ Klingelbergstrasse 82, 4056 Basel, Switzerland 

^ ■ July 4, 2007 

00 ■ 

(N : 

' Abstract 

• The impact of the strongly attractive electromagnetic field of heavy nuclei on electrons in quasi- 

I , elastic (e, e') scattering is often accounted for by the effective momentum approximation. This 

O ' method is a plane wave Born approximation which takes the twofold effect of the attractive nucleus 

^ , on initial and final state electrons into account, namely the modification of the electron momentum 

1 ^ I ' in the vicinity of the nucleus, and the focusing of electrons towards the nuclear region leading to an 

enhancement of the corresponding wave function amplitudes. The focusing effect due to the attractive 

CO ' Coulomb field of a homogeneously charged sphere on a classical ensemble of charged particles 

I incident on the field is calculated in the highly relativistic limit and compared to results obtained 

^ , ■ from exact solutions of the Dirac equation. The result is relevant for the theoretical foundation of 

, the effective momentum approximation and describes the high energy behavior of the amplitude of 

T-H ' continuum Dirac waves in the potential of a homogeneously charged sphere. Our findings indicate 

ly-^ I that the effective momentum approximation is a useful approximation for the calculation of Coulomb 

(•~^ ■ corrections in (e, e') scattering off heavy nuclei for sufficiently high electron energies and momentum 

I transfer 

^p. I Keywords: Coulomb corrections, quasi-elastic electron scattering, effective momentum approxima- 

^> ■ tion. 

, PACS: 11.80.-m Relativistic scattering theory; ll.15.Kc Classical and semiclassical techniques; 

$H ' 1 1.80. Fv Approximations; 1 3. 40. -f Electromagnetic processes and properties; 25.30.Fj Inelastic elec- 

5^ , tron scattering to continuum; 25. 30. Bf Elastic electron scattering; 25.70.Bc Elastic and quasielastic 
scattering 



1 Introduction 

Scattering experiments can be viewed as one of the very important tools of experimental particle physics 
since the famous Lord Ernest Rutherford scattering experiment of a-particles off the nuclei within a gold 
foil in 1911 m. To explore the structure of the nucleus, the main tool used today is electron scattering 
due to the transparency of the nuclear volume for electrons. E.g., inclusive (e, e') scattering, where only 
the final electron is observed, provides information about the nuclear Fermi momentum by measuring 
the width of the quasi-elastic peak 121, or the high-momentum components of nucleon wave functions 
when the tail of the quasi-elastic peak is investigated ||3l|4]|. Information about infinite nuclear matter 
is obtained by extrapolating the mass number A ^ oo f5l, and possible modifications of the nucleon 
form factors inside a nucleus are related to the Coulomb sum rules |6]. However, although electrons with 
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energies of typically some hundred MeV are used in the experiments, the distortion of the electron wave 
functions due to the strongly attractive electrostatic field of heavy nuclei can no longer be neglected, such 
that calculations in the plane wave Born approximation (PWBA) are no longer reliable. 

Calculations using exact Dirac wave functions are feasible but cumbersome and difficult compared 
to the PWBA calculations. As a consequence, various approximate methods have been proposed in the 
past for the treatment of Coulomb distortions |l7l[l|9l[ini[IIl[l2l[l3l[T4l[l5l, and there is an extensive 
literature on the so-called eikonal approximation |n^ [T7l[T8l[T9ll20ll2Tll22ll23l l24 l. 

In this paper, we give a concise classical derivation of the effective momentum approximation 
(EMA), which has the advantage that one works with plane waves and which plays an important role 
in experimental data analysis. The classical high-energy results are compared to results obtained from 
exact solutions of the Dirac equation. Our findings concerning the correct use of the EMA are of actual 
importance, since there is now considerable theoretical and experimental interest in extracting longitu- 
dinal and transverse structure functions as a function of energy loss for fixed three-momentum transfer 
for a range of nuclei. Recently, a Thomas Jefferson National Accelerator Facility (TJNAF) proposal for 
quasi-elastic electron scattering measurement in the momentum transfer range 0.55 GeV/c < |g| < 1.0 
GeV/c was approved such that the experiments will be performed in the near future using ^He, ^^C, ^^Fe 
and ^°^Pb as target nuclei lIBll . 

We shortly comment qualitatively on the connection between the distorted wave Bom approximation 
(DWBA) and the EMA and its correct application. In DWBA, one calculates matrix elements with exact 
initial and final state electron wave functions. Unlike the plane waves with constant amplitude used in 
the PWBA, these wave functions are focused towards the nuclear region, and the local electron momenta 
are enhanced there due to the attractive positively charged nucleus. In the EMA, the focusing and the 
momentum transfer in the relevant nuclear region, where the nucleons get knocked, are accounted for 
by effective (average) values. It is important to mention that one must base EMA calculations indeed on 
average values, although in the literature, the use of effective values for the focusing and the effective 
momenta valid only in the center of the nucleus is widespread. However, the choice of such values is not 
appropriate, as will be explained in detail in this paper. 

There are two equivalent methods for the correct application of the EMA. First, one may calculate 
the cross section from the corresponding theoretical PWBA expression for the (e, e') scattering cross 
section with effective momenta. This introduces an artificially enhanced phase space for the final state 
electron, since also in the DWBA, the phase space is given by the undistorted asymptotic momenta of the 
final state particles. However, this enhanced phase space accidentally accounts for the focusing effect on 
the final state electron with a high level of accuracy. Because the initial focusing has not yet been taken 
into account, one has to multiply the cross section calculated so far additionally by the effective focusing 
factor of the initial state electron. 

Another equivalent approach is to factorize the theoretical expression for the (e, e') cross section into 
the Mott cross section given by eq. (l23l) and a response function according to eq. ((22)) . The interesting 
point is that the impact of the focusing cancels against the modification of the momentum transfer in 
the Mott cross section. Accordingly, one may calculate the EMA cross section by leaving the Mott part 
unchanged and by evaluating the response function with the momenta replaced by their effective values. 

A critical overwiew on the history of the effective momentum approximation and its correct and 
incorrect application can be found in ll26ll . 

2 Quasi-elastic scattering 

In order to illustrate the importance of Coulomb corrections for quasi-elastic (e, e') scattering, we shortly 
review the basic properties of this scattering process. For this purpose, we envisage an electron with 
initial and final asymptotic four-momenta r = {eij,kij), which scatters off a nucleon. We will 
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Figure 1: Quasi-elastic (e, e') scattering cross section data taken at Saclay for initial electron energy 
Ei = 485 MeV and electron scattering angle Bg = 60°. 



always set ?i = c = 1 in the following, and for highly relativistic electrons we have eij = \kij\. 
Additionally, we assume that the nucleon inside the nucleus is at rest. Neglecting interactions on the 
nucleon with its surrounding such that the nucleon can be considered quasi-free, the initial and final 
momenta of the nucleon are given by = {nin, 0) and = {Ef,pf) = (m„ + uj,ki — kf), where 

k^ — k^r is the energy transfer and q = hi — kf the three-momentum transfer of the electron to the 



nucleon. From four-momentum conservation 

we obtain from the four-momentum transfer squared 



2mi 



2mnEf, 



(1) 



(2) 



and consequently lu = [Ef — nin) = Therefore, under the simplifying assumptions made above, 
the (e, e') scattering cross section as a function of the energy transfer for fixed electron scattering angle 
6e should exhibit a peak where 



2mr. 



(3) 



Fig. [T] shows such a typical experimental curve from measurements taken at Saclay ||27]| . First, one 
observes that the peak has a width which is basically due to the Fermi motion of the nucleons. Second, 
the peak is shifted with respect to the empirical formula eq. Q, which predicts Wpeafc — 100 MeV, to 
a value of nearly 140 MeV. A phenomenological description of this observation could be given within 
the Fermi gas model by the observation that eq. ^ does not take into account that an average removal 
energy Ej-em is necessary to remove a nucleon from the nucleus which is larger than the average binding 
energy Ei,ind of a nucleon inside the nucleus. E.g., for ^'^^Pb with Eund — 20 MeV, a two-parameter fit 
for the Fermi momentum kp and the removal energy Erem. leads tokp — 265 MeV and Ej-em — 44 MeV 
|[2l . The higher value of the removal energy also incorporates correlation effects due to the short range 
interaction of the nucleons f2E\. However, there is a significant non-quasielastic background present 
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in Fig. [U which if removed would make the peak appear around 130 MeV, a value which is not so 
different from what one would expect from the binding energy, putting the observations made above into 
perspective. Furthermore, the momentum of the electron in the nuclear vicinity is enhanced due to the 
attraction of the nucleus, which induces an additional positive shift of the peak. This leads us to the idea 
of effective momenta. From a classical point of view, the momentum of a highly relativistic electron 
which moves virtually on a straight line is locally dependent and given by 

= (kij - Vir))kij, (4) 

where kij is the unit vector in direction of kij kij = \kij\, and V{r) is the potential energy of the 
electron in the electrostatic field of the nucleus. This local change of the momentum of, e.g., the incom- 
ing particle with momentum ki = kiki is taken into account by the eikonal approximation through a 
modification of the plane wave part of the free wave function describing the initial state of the particle. 
Defining the relativistic eikonal phase 

z 

Xi{r) = - [ V{r + ks)ds = - [ V{x,y,z')dz' (5) 



if we choose ki = A;* e^, the free electron spinor used in PWBA calculations 

^^{f) = u,^{h)e''^^'' (6) 

is replaced by 

^i{r}=UsMy'^''''^''''^"^ (7) 

in the corresponding eikonal distorted wave Bom approximation (EDWBA). Us^{ki) is the constant 
spinor which depends on the spin (helicity) and momentum of the particle. As desired, the dominant 
longitudinal z-component of the momentum is then recovered via 

The final state wave function is constructed analogously by the replacement e^^f^ — > e^kfT-txf{r) ^ where 



oo 

Xf{^=- j V{r + kfs')ds' . 



(9) 



However, this approximation does not yet include the fact that also the amplitude of the electron wave 
function corresponding to initial and final asymptotic momenta k^j is modified by the attractive nucleus. 
An improved version of the eikonal approximation thus should read 

= /^f (r-)n,,_,(^ij)e*^^-/^^+^>^-/(^"^, (10) 

such that the electron probability density is locally enhanced by focusing factors fij{r). 

A simpler strategy than the eikonal approximation, which will eventually lead to the EMA and which 
avoids the introduction of non-planar wave functions, is to average the locally dependent momentum over 
the nuclear volume, such that effective momenta k'^-^/ are obtained 



Jp{r)d^ 
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with p{r) representing a reasonable nuclear density profile. If both the charge and the nuclear density 
are approximated by a homogeneous distribution inside a sphere with radius R 



I const. : \f\<R 
^(^) = ^ : \f\>R ' ^'^^ 



then it is straightforward to show that the effective momenta are given by 

4, 



^f/ = {hf - ^V{0)) kf = ihf - Veff)kj, (13) 

and the potential energy of an electron 1^(0) in the center of the nucleus is given by ^(0) = 
where a = e^/47r is the fine structure constant and e the elemental charge. Accordingly, one can define 
now an effective four-momentum transfer squared Q^ff- The effective potential V^jf = 4V{0)/5 is 
indeed the average value of the potential V generated by the homogeneous charge distribution inside 
the sphere itself. Eventually, since the present discussion has a phenomenological character due to the 
complex and partially uncertain structure of the nuclear current, one may modify eq. (|3]) to an even more 
general form 

LO = -^+Erem, (14) 

2m„, 



where m„ and Erem are a phenomenological (momentum-dependent) nucleon mass and a phenomeno- 



logical removal energy, respectively. Replacing by Q^ff leads to an additional peak shift of 



MeV in the present example. 



Nucleus 




Figure 2: Quasi-elastic electron scattering off a heavy nucleus. Within a strongly simplified picture, the 
process can be viewed as scattering of the electron off the constituents (mainly nucleons) of the nucleus 
via exchange of a 'hard short-range' photon. 

Viewing quasi-elastic scattering as a nucleon knock-out process provides only a poor picture of the 
actual physical processes taking place inside the nucleus and for details we refer to the literature [29]. 
What is important for the forthcoming section is the fact that the electron interacts with the nuclear 
medium by exchange of photons, and that the hard scattering process can be viewed as a quasi-local 
process. E.g., for ej = 485 MeV and u = 160 MeV, the four-momentum transfer is = (397MeV)^. 
Taking into account that he = 197 MeV fm, the virtuality of the exchanged photon corresponds to a 
typical space-time length scale of 0.5 fm, which is much smaller than the size of the nucleus, as depicted 
in Fig. H 
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3 Effective momentum approximation 



The differential cross section for single nucleon knockout is given by [SO] 

defdnfdEfdQf ~ (27r)' 
with the matrix element 



e}EfPf5{ei+EA-ef-Ef-EA-i)Y^ \Wif\\ (15) 



W,f= d'x d'y / d-'q f^{x) J^{y) 



-iq{x-y) 



(16) 



where Jpf{y) is the nucleon current obtained from some suitable nuclear model, the ^l- G2]l 

indicates the sum (average) over final (initial) polarizations, and Ea, Ea-i is the energy of the initial 
and final nucleus, respectively. In the PWBA, the electron current is given by 

fix) = n,,(^/)7^n,,(^.)e*^^-%n (17) 

where . , Ug^ are initial/final state plane wave electron spinors corresponding to the initial/final electron 

momentum kij and spin Sij. In the DWBA, exact solutions of the Dirac equation are used for electrons 
instead of plane waves. The usual procedure to calculate the inclusive (e, e') cross section is to sum 
over all the individual nucleon knockout cross sections for all protons and neutrons in the nucleus under 
consideration. 

The basic idea of the effective momentum approximation (EMA) is to describe the electron wave 
functions by modified plane waves 

which account for the enhanced electron density and momentum in the nuclear region. Here, k'- ^ and k'l j 
denote effective momenta which need not necessarily be identical. It will be one of the main results of the 
forthcoming section, that for high electron energies and the electrostatic potential of a homogeneously 
charged sphere, k[ j = k" ^ = k^^ is indeed fulfilled, i.e., at high energies, the effective (average) 
focusing factor is given by 

We will therefore identify the k'- j = k'! ^ = k^^ in the sequel. When the exact wave functions appearing 
in the matrix element eq. ( fT6l ) are replaced by the corresponding effective wave functions, the momentum 
integral in eq. (fT6l ) can be trivially performed and replaced basically by a constant factor 1 / q'^ = 
— l/Qgjj. This expresses the fact that the virtual photon emitted by the electron is actually harder than 
if no attractive potential were present, since the electron is accelerated to higher momenta in the nuclear 
vicinity, and ^ is the photon virtuality averaged over the nuclear volume. 

Note that the reason why the replacement of the locally dependent wave function amplitudes and mo- 
menta by effective values makes sense is rooted in the local character of the scattering process mentioned 
above. E.g., if the virtual photon would propagate over distances comparable to the size of the nucleus, 
then nucleons could also be knocked by photons which were emitted outside the nucleus, such that an 
averaging of the focusing and the local momenta inside the nuclear interior would not make sense. The 
correct mathematical counterpart of this pictorial description can be found in HI. For the EMA to hold, 
it is mandatory that the wave lengths of the electron and the virtual photon are significantly smaller than 
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the nuclear radius, i.e. e/ > 200 MeV and > (200 MeV)^ corresponding to a length scale of 1 fm 
should be required for ^osp^, |[3Tll32l . Note also that the enhancement of the wave function amplitudes 
is not very large at high energies. One can write 

fl//ir) = 1 + Sijir), \5{f)\ « 1, (19) 
and we may therefore neglect higher order terms in the 6's in formal expressions, like 

(///^(f)/;/^(r)) = ((1 + 6ij{fm + kfim ^ {fl^\mfy\r)), ao) 

and one may equate expressions like 

fijif) = (1 + 6ij{r}f ~ 1 + 26ij{r). (21) 

It is instructive to calculate the impact of the focusing factors on the size of the cross section for a 
typical example. E.g., if we consider electron scattering off ^'^^Pb for \ki\ = 485 MeV/c and lo = 100 
MeV, we have V{0) = -25 MeV and 4F(0) /5 = -20 MeV such that kf^ is given by (485 + 20) MeV, 
and k^-^ = (385 + 20) MeV. The focusing factors enter the cross section both linearly via the matrix 

element squared, enhancing the cross section by a factor of {k^^^ /ki)'^{ky^ /kj)"^ = 1.2. 

There are two different, but equivalent strategies to calculate cross sections in the EMA framework. 
First, the EMA cross section can be calculated by replacing the electron momenta kij by the effective 
momenta k'^^ in the (theoretical) expression for the quasi-elastic scattering cross section (accordingly, 
the energies e^j must be replaced by \k^y\)- The cross section obtained this way must be multiplied 
subsequently by the factor (k'^/ki)"^ which accounts for the focusing of the incoming electron wave in 
the nuclear center. The focusing factor (k'j/kf)'^ for the scattered electron is already contained in the 
artificially enhanced phase space factor of the final state electron, if k'j = k'j is presumed. Second, the 
cross section for inclusive quasi-elastic electron scattering can also be written by the help of the total 
response function Stot as 

^-^f^ = <TMott X Stot{\q\,i^,&e), (22) 

where the Mott cross section is given by (g^ = Q^) 

TMott = cos\Qe/2)e} /qf,. (23) 

The Mott cross section remains unchanged when it gets multiplied by the EMA focusing factors and the 
momentum transfer is replaced by its corresponding effective value. A short calculation shows indeed 
that (eij » m) 

^eff ^ ^ ^/ J^^^f ) ^ ^ ^24) 

kikf ^eff 

Therefore, the EMA cross section can also be obtained from (|22]) by leaving the Mott cross section 
unchanged and by replacing Stot{\q\,^, ©e) by the effective value 

Stot{\qeff\,^,Qe) = Stot{\kl^^ - ky^\,LL>,Qe)-, (25) 

since the effect of replacing q^^ by its effective value in the Mott cross section is to exactly divide away 
the initial state focusing factor and the final state focusing factor which is generated by the replacement 

e// l;e//| 

e J — > e J = \kj \. 
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Nucleus 



Figure 3: Electrons incident on an nucleus with impact parameter 60 ■ 



4 The classical focusing factor 

The focusing factor can be derived approximately from a classical toy model according to Figs. |3]and 
m We consider the trajectories of an ensemble of highly relativistic particles approaching a nucleus 
located in the center of the (6, z)-coordinate system. The particles shall move in z-direction with equal 
velocity and with an asymptotic impact parameter in the range between 60 and 60 + dbQ. The longitudinal 
velocity of the particles can be taken as the speed of light, since the particles are highly relativistic and 
changes of the velocity in transverse direction and the kinetic energy cause a negligible second-order 
effect to the longitudinal component. Therefore, one may adopt the straight-line approximation by setting 
z{t) = ct = t, r{t) = y^Q + t^, and the impact parameter will be considered constant at certain stages 
of our calculation. Due to the attractive nucleus, the original impact parameter 69 is reduced to 6(60 ) z) 
as a particle moves along its trajectory, such that the particle density at z is increased by a focusing factor 
/ which is given by the ratio of the area of two annuli with radii ^0 + dh^ and 6(60, z),h{hQ + dbo, z): 

f {bo,z) ~ / {b{bQ,z),z) = — — . (26) 

OOq Oq 

In the following, we will calculate b{bo,z) for an electron in the potential of a homogeneously 
charged sphere with radius R and charge eZ, given by 

I : r>R 

since the potential of a homogeneously charged sphere provides a simple but quite realistic model for the 
electromagnetic field of a heavy nucleus like, e.g., ^'^^Pb, where one has i? ~ 7.1 fm and Z = 82. 
The force ft acting on the particle in transverse direction is given by 

r bdVhom{r) bodVhomir) 

ft = 5 ' ^ , (28) 

r or r or 

In eqns. ( [261 ) and ( |28] ). we made use of the straight-line assumption by replacing b by b^. Also in the 
forthcoming, we will sometimes replace b by 69 or use these two quantities synonymously where such 
a substitution is adequate. To calculate the transverse acceleration of the particle due to the attractive 



8 



I 




I 

Figure 4: The three different regions according to the case distinction in the text. 



Coulomb field, we distinct three cases (see Fig. H)). In the first case, we consider the region where the 
particle moves solely in the 1/r-field according to the straight-line approximation, i.e. where we have 
bo > R (or 60 < —R^ if we formally allow negative impact parameters), plus the region with bo < R 
and z < — — 6q. The transverse force is then given by 



aZbo 



(29) 



Correspondingly, we obtain for the transverse acceleration, taking into account that the "transverse mass" 
of a relativistic particle is given by its energy E, which is also considered as constant: 

aZ bo aZ 69 



where r is the distance of the particle from the nuclear center, and 



(30) 



vi{bo,t)= / al{bo,t')dt' 



aZ t + yjt"^ + i>o _ aZt + r 
boV^T^ ~ ~^~b^' 



(31) 



Note that from vl{t ^ 00) = — we obtain for a pure 1/r-Coulomb field the well-known transverse 
momentum transfer 

= (32) 



ho 

Furthermore, we obtain from ( [3T]) 

z 

b\bo,z) = bo+ / vf{bo,t')dt' = bo 



aZ z + r 



(33) 



A short calculation yields the focusing factor 



aZ 11 z + r 



1 -1 



E 



+ 




-1 



(34) 



Since we are interested in the high energy behavior, we keep in eq. (l34l ) only the relevant zeroth- and 
first-order terms in aZ/E. For the focusing factor in region I we obtain the simple result 



/ (60,^) = 1 + —-. 



(35) 
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The calculation for regions II and III are a bit more involved, but can be performed along the same lines 
as above. 

We calculate now the focusing inside the charged sphere (region II), which gets traversed by particles 
with bo < R. Inside the sphere, the transverse acceleration of the particles is due to the harmonic 
oscillator potential generated by the homogeneous charge distribution. Correspondingly, we have 

haZ r aZ bo 
=—r^R-^^-^R-^' ^^^^ 
and for the transverse distance from the z-axis we obtain after a short calculation 

~N 2 



b''{bo,z) = bo-^ 



R-R ^ 'f[^ + R) ^ [R-R}[R + z 



bo R^ boR 



(37) 



Above, we have introduced the abbreviation R = ^R? — 5q. Note that the first term in the bracket 
above describes the transverse shift of the particles when they arrive on the surface of the charged sphere 
according to eq. (l33l) . where z = —R. The second term is due to the transverse acceleration of the 
particles inside the sphere, and the last term in the expression above is generated by the transverse 
velocity vj:{bo,—R) which is reached by the particles when they cross the border of the sphere. A 
straightforward calculation leads to the following result for the focusing factor at first order in aZ/E: 

f\ho, = 1 + + 3^ + ^ - 2bl/R') . (38) 

In the center of the nucleus, the particle density is enhanced by a factor 

= f\bo = 0, z = 0) = 1 + (39) 

however, the average focusing factor inside the sphere is given by the volume integral 

/ f'{bo.z)dV 

i.e. one obtains the focusing factor used in the effective momentum approximation, where the increased 
particle probability density near the nucleus is taken into account by multiplying the particle's Dirac 
wave function by a suitable factor /^/^ ~ {P^Y^'^- 

Finally, we consider the 'shadow region' of the nucleus (region III). We calculate first the transverse 
velocity of the particle when it arrives in region III in three steps. First, the particle moves inside the 
1 / r-field and reaches a transverse velocity 

T ~, aZ R — R 
v,=vi{bo,-R) = - — -^ (41) 

at the surface where it enters the sphere. Inside the sphere, the particle undergoes a constant transverse 
acceleration a^^ for t G [—R, R]. Therefore, the particle gains an additional velocity 

aZ2boR 

V2 = (42) 

In the downstream region III, the particle is again moving in the field of a point-like charge, and the 
transverse acceleration is given by 

ar = -^^^, (43) 
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such that we end up with (r = \/{t^ + ^o)) 



V1 + V2 + 



j al"{bo,t')dt' 



R 



aZ 



R 1 26oi? 



aZ 

~e' 



t + r 



R 



aZ 



t + r 2R 



+ 



bor boR 
2boR 



1 

b'o 



bor boR ' R^ 

The transverse distance of the particles from the z-axis in region III is therefore given by 



(44) 



b''\bo,z)=b'\bo,R)+ I vi''{bo,t)dt 



.111/ 



R 



aZ 



'bo R bo R^ 



+ I vl'\bo,t)dt 



R 



aZ 



z 2bozR j_ 52 2zR 

Vo ^ R^ ^ 60 kiR 



For the focusing factor we obtain from eq. 



f'\bo,z) = l+'^ 



6z 



6zbl 



RR y/z^ + b'^ R^R 



(45) 



(46) 



We finally summarize the results as follows. An attractive nucleus modeled by a homogeneously 
charged sphere acts like a focusing lens on an ensemble of classical particles incident on the nucleus with 
impact parameter 60 on quasi-straight trajectories parallel to the z-axis. For highly relativistic particles, 
the particle density is enhanced by a focusing factor 



fibo,z) = l + ^^bo,z) 



(47) 



with (r = ^/Jbl + z^)) 



$(60, z) = < 



3 I Q2-R I 



2bl/R^ 



1 I 6z 

r ^ RR R3R 



I 

II 
III 



(48) 



The typical deviation of the focusing from unity in the nuclear interior is of the order of 3aZ/ER or 
Vhom{0) / E; this ratio should be considered as the expansion parameter for higher order corrections to 
the focusing, which become irrelevant at high energies. 

Fig. [5] shows a surface plot of the universal function $(601 z). In order to compare the classical 
focusing to the results obtained by solving the Dirac equation exactly ll32l[33l[34l . we define 



HE;b,z) :=^{p{E;b,z)-l) 



(49) 
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Figure 5: Left: The function $ describing the focusing effect of the attractive potential of a homoge- 
neously charged sphere on an ensemble of highly relativistic particles. In order to symmetrize the figure, 
the impact parameter b can also be negative. Parameter values typical for a '^^^ Pb nucleus have been used 
{R = 7.1 fm, Z = 82). Right: Focusing function ^{E; b, z) for an electron incident on the potential with 
E = 200 MeV and positive helicity. The left and the right plot agree well inside the sphere, however, 
the wave focusing is clearly larger on the rear side of the nucleus than in the classical highly relativistic 
case. 



where p{E;b,z) = "i^{E;r)^^{E;r) is the axially symmetric probability density of the Dirac wave 
function "^{E; f) of an electron incident with asymptotic momentum k = \/ E'^ — m?z ~ Ez and spin 
parallel to the z-axis. Results are shown in Figs. [5] and |7j The Dirac density indeed approaches the 
classical limit for high electron energies. Note that the focusing of the Dirac wave function is clearly 
underestimated by the classical high-energy approximation in region III, where the straight-line assump- 
tion starts to break down. However, in the case of quasi-elastic electron scattering, the relevant region 
is the interior of the nucleus, where the focusing is described in a satisfactory way for electron energies 
above 200 MeV as shown in Figs. |5l[6]and|7] 

We note that Knoll [8 | derived the focusing effect from a high energy partial wave expansion, fol- 
lowing previous results given by Lenz and Rosenfelder CKTOl. For the incoming particle wave expanded 
around the center of the nucleus he obtained 

{1 + air^ - 2a2kif+ iair'^kiV + ia2[{ki x rf + a{ki x f)]}usi{ki), (50) 

where 5i is a phase, ki is an effective momentum parallel to ki calculated by using the central potential 
value ki = ki — Vhom{^)^ and a acts on the spinor Us^ih) to describe spin dependent effects, which are 
negligible in our cases of interest with definite helicity. An analogous equation holds for the distortion 
of the outgoing wave. The parameters 01^2 depend on the shape of the potential. For a homogeneously 
charged sphere with radius R they are given by 

aZ 3aZ 

ai = = , 02 = ^ , (51) 

6kiR^ 4k? R^ 
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-5 -5 
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Figure 6: The ratio ^{E, b, z)/^{b, z) for E = 200 MeV in region n. 




Figure 7: Left: Longitudinal focusing of the classical particle density $(6 = 0, z) versus the focusing of 

the Dirac wave function ^{E, b = 0, z) for two different electron energies E = 200 MeV and E = 400 
MeV. Right: Corresponding transverse focusing ^{b, z = 0). The bar at the bottom depicts the extension 
of the nucleus. 
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and the central potential value is given by Vhom{0) = ~f The increase of the amplitude of the wave 

while passing through the nucleus is described mainly by the — 2a2A;ir-term, the oir^-term accounts 
for a decrease of the focusing also in transverse direction. Performing the replacements = 6^ + z^, 

kif = kiZ and | A;^ x f| = kib, one obtains from eq. (l50l ) for a particle with spin parallel to the momentum 
and energy E = ki up to second order in b and z and first order in aZ 



*i(f)T^i(f 



aZ 



Us^{ki)^Us^{ki) 



(52) 



i.e. the correct linear term in z in eq. (|48T | is recovered, however, the transverse decay of the focusing 
is strongly suppressed in the V^ff = 4y(0)/5 expansion given by Knoll, which is therefore not suited 
to describe the focusing for fe 7^ 0. The use of the expansion eq. (l50l ) is the reason why the (e, e') cross 
sections in [23] are overestimated. 

Calculations using exact Dirac wave function show that the effective electron momenta are very well 
described by an effective potential V^ff = 41/(0)/5 fST]. Furthermore, our findings indicate that also the 
focusing can be described in an accurate way by the same effective potential value. This demonstrates 
the validity of the EMA as a valuable tool for the description of Coulomb distortion effects. One may 
observe from Fig. |7]that the exact focusing for finite energies is slightly larger than in the high energy 
limit. This leads to a minor amplification of the DWBA cross section compared to the EMA result. 
Exact calculations, which will be presented in a forthcoming paper, show that this effect is only of the 
order of 2% in the region of the quasielastic peak for typical kinematical values used in experiments, e.g. 
for an initial electron energy of = 485 MeV and scattering angle 0e = 60°, or = 310 MeV and 
Ge = 143°. 



5 Final remarks and conclusions 

The high energy trajectory of a charged classical particle moving in the field of a homogeneously charged 
field was investigated and related through the quantum-classical correspondence principle to the proba- 
bility density of exact continuum wave functions obtained as solutions of the Dirac equation. As a result, 
a universal function $ was found which allows to describe the high energy behavior of the amplitude of 
Dirac electron wave functions with definite helicity. The focusing in the downstream side of the charge 
distribution converges slowly towards the high-energy limit described by $, however, the universal func- 
tion <I> provides an accurate description of the focusing inside the charged space region, which can be 
considered as a model for the charge distribution of a heavy nucleus. As a consequence, it is found that 
both the effective (average) momenta and the average focusing can be described by a common effective 
potential Veff = 4y(0)/5 in the case of a homogeneously charged sphere, despite the fact that the local 
classical momenta exhibit the same spherical symmetry as the electrostatic potential, whereas the axially 
symmetric focusing is smaller in the upstream side and larger in the downstream side of the nucleus. 

Our findings establish the role of the EMA as a valuable semiclassical method for the analysis of 
Coulomb corrections in (e, e') scattering. They also indicate that the analysis of experimental data based 
on calculations of Kim et al. should be revisited |[T5l l35l . and support the strategy in previous works 
concerning the extraction of the longitudinal and transverse response functions in medium-weight and 
heavy nuclei |^]. However, it is also advisable to await new experimental data which will hopefully be 
accessible in the near future ll25l . 

We finally remark that Baker investigated also the second-order eikonal approximation for potential 
scattering in the non-relativistic case [36], finding thereby an expression for the focusing factor of con- 
tinuum Schrodinger wave functions. For the focusing in the center of a spherically symmetric potential. 
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one finds (see eq. (23) in 1361) 



/V2(o) = /i/2(6 = o,. = 0)=.l-^, or /(O):.!-^ (53) 

Zkv kv 

wliere k is tlie asymptotic momentum and v the velocity of the particle. Roughly speaking, the approx- 
imation is valid if the kinetic energy of the particle is larger than the depth of the disturbing potential 
m S> -Efcin = E — m ^ ^(0), and the wave length of the particle ~ 2iT/k should be significantly 
smaller than the extension of the potential. For the classical particle momentum in the center of the 
potential fc(0) one has non-relativistically 



k{0) = V2m(i?fc,„ - y(0)) = ^/2mEk,n\\ 1 - ^ =^ A: 1 - , (54) 

such that 

,55) 

kv k 

i.e. it is found that the probability density is enhanced by the ratio of the central and asymptotic momenta 
k{0)/k, instead of {k{0)/k)^ ~ {{E - V{0))/Ef in the highly relativistic case. One may ask how the 
non-relativistic and the highly relativistic regime are connected. A classical relativistic analysis of the 
particle trajectories shows that the central focusing is given by the expression 

;(„).M^. ,5« 

which interpolates between the non-relativistic and relativistic regime and which is given here, for the 
sake of brevity, as a result without proof. 
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